Observation of Anomalous Spin Segregation in a Trapped Fermi Gas 
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We report the observation of spin segregation, i.e., separation of spin density profiles, in a trapped 
ultracold Fermi gas of ^Li with a magnetically tunable scattering length close to zero. The roles of 
the spin components are inverted when the sign of the scattering length is reversed. The observed 
density profiles are in qualitative agreement with the spin- wave theory applied previously to explain 
spin segregation in a Bose gas, but disagree in amplitude by two orders of magnitude. The observed 
atomic density profiles are far from equilibrium, yet they persist for ^ 5 seconds in a trap with 
an axial frequency of ^ 150 Hz. We attribute this long lifetime to Fermi statistics: The scattering 
amplitude is nonzero only for atoms in opposite states, and vanishes for atoms in the same state. 
By measuring the magnetic field at which spin segregation ceases, we precisely determine the zero 
crossing in the scattering length of ^Li as 527.5 it 0.2 G. 



Ultracold Fermi gases near a Feshbach resonance [1] ex- 
hibit strong interactions [2] and offer unprecedented op- 
portunities to test nonperturbative quantum many-body 
theory in systems from high temperature superconduc- 
tors to nuclear matter. Near the resonance, a bias mag- 
netic field tunes the interaction between two spin states 
from strongly repulsive to strongly attractive. These 
strongly interacting Fermi gases have been extensively 
studied over the last few years [3]. 

In contrast, the weakly interacting regime of a Fermi 
gas near a Feshbach resonance has received relatively lit- 
tle attention. For a bias magnetic field near the zero 
crossing, the s-wave scattering length for atoms in oppo- 
site spin states can be tuned smoothly from small and 
positive to small and negative as the bias field is varied. 

We report the observation of spin segregation in an 
optically-trapped ultracold Fermi gas of ^Li with a bias 
magnetic field close to the zero crossing. A sample of ul- 
tracold ^Li gas is created in a coherent superposition of 
the two lowest spin states |1) and |2), with equal popula- 
tion in each state. Initially, the two spin states have the 
same atomic density profile. After several hundred ms, 
the spatial densities segregate, with one spin state mov- 
ing outward and one moving inward, as shown in Fig. [H 
The roles of the spin states are interchanged when the 
sign of the scattering length is reversed. The differential 
trapping potential experienced by atoms in the two states 
cannot account for the spin segregation, since it is four or- 
ders of magnitude smaller than the thermal energy ksT. 
A longitudinal spin wave theory 0, [HI, 0] has been used to 
explain spin segregation that was observed in a trapped 
Bose gas [7]. In the Bose gas, numerical simulations are 
in good quantitative agreement with the experimental re- 
sults. In contrast, our observations of spin segregation in 
a Fermi gas show great discrepancy with the theoretical 
calculations based on a spin wave assumption. 

In the experiments, a sample of ^Li atoms in a 50/50 
mixture of the two lowest hyperfine states is loaded into 
a CO2 laser trap with a bias magnetic field of 840 G, 
where the two states are strongly interacting. Evapora- 
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FIG. 1: Absorption images (state 1 and state 2) taken at 
200 ms after the RF pulse for 526.2 G (scattering length a < 
0) and for 528.8 G (a > 0). Each image is 1.2 mm in the 
horizontal direction. 



tive cooling is performed to lower the temperature of the 
sample [2]. The magnetic field is then increased in 0.8 
seconds to a weakly interacting regime at 1200 G where 
an on-resonance optical pulse of 15 /is is applied to re- 
move atoms of one state while leaving atoms in the other 
state intact. With the single state present, the magnetic 
field is lowered in 0.8 seconds to 528 G, near the zero 
crossing. Then an RF pulse is applied on the |1) — |2) 
transition to create spin coherence. Instead of applying 
a 7r/2 pulse, the RF frequency is swept across the reso- 
nance of the hyperfine transition ~75.60 MHz to achieve 
the coherent spin transfer. We control the amplitude of 
the spin transfer by varying the sweep rate. Typically, we 
sweep the RF frequency by 35 kHz (centered at hyperfine 
transition) in 40 ms in order to transfer about 50% of the 
atoms. Finally, we take absorption images of atoms in 
both states (in separate experimental cycles) at various 
times after the RF pulse. From those images, we can de- 
termine the parameters of the cloud sample such as tem- 
perature and density. At the final optical trap depth, the 
measured trap oscillation frequency in the transverse di- 
rections is ij± = 27rx4360 Hz, while the axial frequency is 
iJz = 27r X 145 Hz at 528 G. Atoms in each of the two spin 
states experience slightly different magnetic trap poten- 
tials in the bias magnetic field. Due to the bias magnetic 
field curvature, the differential potential Ui — U2 varies 
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across the sample and (P{Ui - U2)/h/dz'^ = 4.4 x 10"^ 
Hz//im^. The total number of atoms is :^ 2.0 x 10^. 
The corresponding Fermi temperature is :^ 7 jaK. 

We find that the spin segregation is not sensitive to the 
temperature of the cloud. We observed spin segregation 
over the range from a degenerate Fermi gas (T O.lSTp) 
to a non- degenerate one (T 4Tf). All the data shown 
in this Letter were taken in the non-degenerate regime. 

The experimental results for spin segregation as a func- 
tion of bias magnetic field are summarized in Fig. [2l 
Here, the sample temperature is T 27 fiK and the 
peak atomic density is 1.2 x 10^^ /cm^. The axial 1/e 
width for a fit of a gaussian distribution to the initial 
density profile of the sample is 300 /im. 

We use the data at 529.3 G, Fig. [21 as an example 
to describe the observed spin segregation in our experi- 
ments. Note that magnetic fields are calibrated with an 
RF spectroscopic technique, as described below. Spin 
segregation occurs only after a coherent superposition 
of the two states is created. Spin segregation starts to 
build up during the RF pulse duration of 40 ms, as we 
sweep the RF frequency across the hyperfine transition. 
Starting about 40 ms after the RF pulse, we observe a 
significant change in the spatial profiles of the initially 
overlapping clouds for states |1) and |2). The cloud size 
of state |1) decreases while the density profile of state |2) 
evolves into a double-peak structure. Maximum segrega- 
tion occurs at about 200 ms and is maintained for several 
seconds. Then the two segregated clouds relax back to 
their original overlapping profiles in ~9 seconds. 

Spin segregation observed in an ultracold Bose gas 
has been described previously as overdamped spin 
wave 0, H, Q . The basic idea can be described in terms 
of a Bloch vector. We take the axis w as the "longitudi- 
nal" population inversion, and the axes u and v as the 
"transverse" coherence. An RF pulse prepares a coherent 
superposition of spin-up and spin-down with polarization 
in the u-v plane. Due to the magnetic field inhomogene- 
ity, the hyperfine transition frequency varies along the 
long axial direction of the sample ~10 Hz. After the RF 
pulse, the inhomogeneous precession of the Bloch vectors 
about the w-axis starts to build up a spin orientation gra- 
dient in u-v plane along the axial direction. Meanwhile, 
atoms with different spins move around and collide with 
each other. The interaction due to the binary collision 
leads to precession of each atom's spin about the total 
spin vector of both atoms. Since both atoms have spins 
in u-v plane, the precession of each spin about the to- 
tal spin in u-v plane produces spin components out of 
the u-v plane. The subsequent formation of di w com- 
ponent leads to a spatially varying population inversion 
and accounts for the spin segregation. 

We tried to use a spin-wave theory to explain the 
spin segregation observed in our experiment. We derived 
equations for the spin density and spin current using a 
Wigner function operator approach fid\. The Wigner 



phase space operators are written in terms of position 
space creation and annihilation operators 

1) ■(') 

The Heisenberg equations are Waa' (r, p) = 
|- [ff, VFq,c^/ (r, p)], where H is the Hamiltonian op- 
erator 
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j(fvijd^Y2^{vi)ip+{T2)Uju{vi-T2)'lpk{y^2 



with Uj — \vfHjj\{x^ +^^) + arid '^j7c(ri — = 

^figjkS{ri — r2). Here, j, k refer to states |1) or |2); Uj is 
the total external potential (combined optical and mag- 
netic) for atoms in state j; hcoj is hyperfine energy of 
state j; gjk = 4:7rhajk/m with m thes atomic mass, ajk 
is the scattering length for atoms between state j and 
state k. 

Spin density operators, Sx = [ipi {r)ip2{'r) + H.c.]/2, 
etc., can be written in terms of the momentum integral 
of the Wigner phase space operators. 
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The corresponding spin current operators are obtained 
using 

ih 

m 



We obtain the spin- wave equations for both a Bose gas 
and a Fermi gas. 



dt 



(4) 
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vs^'^ = {n' xj)^'^ -7J 



2m 4m ' ^ ^ 



where 



5uj2 = 2giini - 2^22^2 - 2(e + l)gi2S^^\ 
= n^2egi2S, 

U = [/i + [/2 + 2%iini + 2%22^2 + %12^, 

= hn\ 

Here i = x^y^z is the spin space index: S*^*^ is a spin 
density. For example, S'^^^ = h{ni — 712) with ni, n2 
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FIG. 2: Axial atomic density profiles (integrated over transverse directions) at various times after the RF pulse for different 
magnetic fields. The scattering length is estimated with a{B) = d{B — Bo), where d = 3.5 ao/G [1, 0| and Bo = 527.5 G 
(described below) Note that different graphs have different time scales. 



the atom densities in states |1) and |2). J*^*^ is the spin 
current corresponding to spin component i [11]. ujo/^tt is 
the hyperfine transition frequency, a = ksT /m. nis the 
total atom number. 7 is the relaxation rate, which we 
take to be in modelling the data. Finally, e = +1(— 1) 
for Bosons (Fermions) and an = a22 = for Fermions. 

Numerical simulations based on these equations quan- 
titatively reproduce prior calculations of spin segregation 
in Bose gases [2, 'H', ^ that agree with the experiment. 
In contrast, these equations only qualitatively describe 
the spin segregation observed in our experiments with 
a Fermi gas. The equations correctly predict the inver- 
sion of the roles of the spins when the scattering length 
is reversed in sign. However, the predicted amplitude 
of the spin segregation derived from the calculations is 
two orders of magnitude smaller than what we observe; 
the predicted segregation time is one order of magnitude 
smaller than what we observe. This suggests that addi- 
tional physics is important for explaining spin segregation 
in a Fermi gas. 

Fig. [2] also shows the damping time (defined as the 
time for the cloud to relax back to its original density 
profile) increases as the scattering length is reduced in 
magnitude. It is interesting to note that for certain mag- 
netic fields (for example, 526.2 G, where the scattering 



length is about 5ao, ao is Bohr radius), spin segregation 
can be maintained for up to 5 s. This result is strik- 
ingly longer than that observed in the previous experi- 
ment with a Bose gas, where the samples relaxed from 
maximum spin segregation back to their original density 
profiles in about 200 ms. We attribute this result to the 
facts that ultracold fermionic ^Li atoms interact only via 
s-wave scattering between opposite spin states and that 
the ^Li scattering length near the zero crossing is much 
smaller than for the Bose gas experiments, ~ 5ao com- 
pared to - lOOao for ^^Rb. 

We also note that the decoherence time in our experi- 
ment is measured to be about 70 ms, which is compara- 
ble to the longest spin segregation time that we observe. 
The decoherence time is the measured time scale for the 
phase of the coherent superposition to decay, producing 
a simple binary mixture. Decoherence is caused by the 
spatial variation of the magnetic field across the sample. 
The short decoherence time may explain why we can not 
clearly observe spin segregation when the magnetic field 
is too close to the zero crossing. When the scattering 
length is very small, the time required for spin segrega- 
tion to build up is longer than the decoherence time. 

An additional experiment was performed to further 
study the spin segregation. After we create a sample 
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FIG. 3: Axial atomic density profiles at 200 ms after coher- 
ence is created by an RF pulse, for different magnetic fields. 
The zero crossing point in the scattering length is determined 
by studying the [transition-NOT a transition] in the spin seg- 
regation behavior as the magnetic field is varied. 

that reaches the maximum spin segregation in 200 ms at 
529 G, we immediately blow out atoms of one state using 
a resonance optical pulse. Then we take images of the 
remaining state at various times after the optical pulse. 
We observe that spin segregation maintains its atomic 
density profile for ~ 2 seconds. This shows that slow re- 
laxation of one state does not depend on presence of the 
other state. As we know, a single-peak gaussian distribu- 
tion is a stationary equilibrium state for a thermal cloud 
in a harmonic trap. In contrast, the double-peak and 
one-peak non-gaussian atomic density profiles observed 
in our experiments are far from equilibrium. It is re- 
markable that we can create a non-equilibrium stationary 
system that has such a long lifetime, ~ 2 s, compared to 
the axial trap period 6.9 ms and segregation time ~ 200 
ms. 

Previous experiments have shown that a zero crossing 
occurs in the |1) — |2) scattering length of ^Li near 528 
G 0, [Til . We now describe a new determination of the 
zero crossing based on the behavior of the spin segrega- 
tion as the magnetic field is varied. As described above, 
the roles of the spin states are interchanged as the scatter- 
ing length a changes sign from negative to positive. Fig. [2] 
shows data for a < at 526.2 G and for a > at 529.3 
G. The spin segregation time (defined as the time for 
the cloud to reach maximum spin segregation) decreases 
(increases) as the scattering length increases (decreases). 
Spin segregation becomes difficult to observe when the 
magnetic field is too close to the zero crossing, for exam- 
ple, 527.3 G in the figure, as the segregation time then far 
exceeds the decoherence time, as noted above. Similarly, 
very far from the zero crossing, ~15 G above or below 
collisional relaxation dominates and spin segregation is 
not observable (not shown in the figure). 

Taking advantage of the fact that spin segregation 
slows for magnetic fields very close to the zero crossing, 
we determine the zero crossing precisely. We take ab- 
sorption images of both states at 200 ms after the RF 
pulse for different magnetic fields. Fig. [3] shows a clear 
change in the spin segregation behavior from 526.2 G to 
529.3 G. The magnetic fields are calibrated by RF spec- 
troscopy using the Hamiltonian for the ^Li ground state 
hyperfine interactions in a magnetic field, 

H,nt = ^S-l- ^{gjS + ^,1) . B, (6) 



where /j^s/h = 1.399624604 MHz is the Bohr magneton, 
ahf/h = 152.1368407 MHz [14] is the ground state mag- 
netic hyperfine constant for ^Li, gj = —2.0023010 jlj] 
is the total electronic g-factor for the ^Li ground state, 
gi = 0.0004476540 [14] is the nuclear g-factor, h is 
Planck's constant, and B is the bias magnetic field. 

The transition in the behavior of the segregation occurs 
at 527.5 G, which is determined as follows. We fit the 
axial atomic density profiles of both states with a gaus- 
sian distribution for the data taken at 527.3 G and 527.8 
Gs. For each data set, we derive gaussian widths of both 
clouds and their ratio. The ratio of the width of state 
|1) to that of state |2) is greater than 1 for 527.3 G; less 
than 1 for 527.8 G. By interpolation, we find the point 
where the ratio of widths is equal to 1, which determines 
the value of magnetic field 527.5 ±0.2 G (the error bar is 
the quadratic combination of the shot-to-shot fluctuation 
of the magnetic field setting and the magnetic field cal- 
ibration error arising from the RF spectrum fit). Since 
no spin segregation occurs with the scattering length of 
zero, we interpret the value of magnetic field where the 
ratio of widths is 1 to be the zero crossing. 

In conclusion, we have observed spin segregation in a 
trapped ultracold atomic Fermi gas of ^Li with a scat- 
tering length close to zero. The observed behavior of the 
spin segregation appears to require a modification of the 
current theory based on spin waves, or possibly a new 
mechanism. In addition, using spin segregation, we have 
precisely determined the zero crossing in the scattering 
length of ^Li. 

This research is supported by ARO, NSF, and DOE. 



* jet@phy.duke.edu 

[1] E. Tiesinga, A. J. Moerdijk, B. J. Verhaar, and H. Stoof, 
Phys. Rev. A 46, R1167 (1992); E. Tiesinga, B. J. Ver- 
haar, and H. Stoof, Phys. Rev. A 47, 4114 (1993). 

[2] K. M. O'Hara et aL, Science 298, 2179 (2002). 

[3] S. Giorgini, L. P. Pitaevskii, and S. Stringari, arXiv: 
0706.3360 (2007), to appear in Rev. Mod. Phys. 

[4] M. 6 Oktel and L. S. Levitov, Phys. Rev. Lett. 88, 
230403 (2002). 

[5] J. E. Williams, T. Nikuni, and Charles W. Clark, Phys. 

Rev. Lett. 88, 230404 (2002). 
[6] J. N. Fuchs, D. M. Gangardt, and F. Laloe, Phys. Rev. 

Lett. 88, 230405 (2002). 
[7] H. J. Lewandowski et a/., Phys. Rev. Lett. 88, 070403 

(2002). 

[8] M. Bartenstein et al, Phys. Rev. Lett. 94, 103201 (2005). 
[9] P. S. Julienne, private communication. 
[10] J. E. Thomas and L. J. Wang, Phys. Rep. 262, 311 
(1995). 

[11] J is a tensor and J Ms a vector in coordinate space. 
[12] K. M. O'Hara et a/., Phys. Rev. A 66, 041401(R) (2002). 
[13] S. Jochim et al, Phys. Rev. Lett. 89, 273202 (2002). 
[14] E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. 
Phys. 49, 31 (1977). 



